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ABSTRACT 

The  factors  affecting  the  stability  and  the  maneuvera- 
bility of  a  dynamically  unstable  ship  were  studied  using 
the  linear  and  non-linear  equations  of  the  motion  of  the 
ship.   The  DSL/360  Language  was  used  to  simulate  the  dynam- 
ics and  to  study  both  unstable  motions  and  standard 
maneuvers. 
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Note 


1)  Symbol  of  prime  ( ' ) ,  except  where  otherwise  stated 
indicates  the  values  of  an  item  non-dimensionalized. 

2)  Symbol  of  dot  (•)  indicates  differential  value  of 
an  item  with  respect  to  time. 

3)  Subscript  f  indicates  the  value  of  an  item  of  fin. 

4)  Subscript  h  indicates  the  value  of  an  item  of  ship's 
hull. 

5)  Subscript  f,  indicates  the  value  of  an  item  of  the 


large  deadwood  included  rudder. 

Subscript  f?  indicates  the  vai 
small  deadwood  or  rudder  only 


6)   Subscript  f~  indicates  the  value  of  an  item  of  the 
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I.   INTRODUCTION 

A  ship  is  said  to  be  dynamically  stable  on  a  straight 
course  or  in  a  turn  of  constant  curvature  if,  when  slightly 
disturbed  from  it's  steady  motion,  it  will  soon  resume  that 
same  motion  along  a  slightly  shifted  path  without  any  cor- 
recting control  being  applied.   A  ship  which  is  dynamically 
unstable  in  straight  line  motion  cannot  maintain  straight 
line  motion  when  there  is  no  control.   However,  with  regard 
to  maneuverability,  a  ship  which  is  too  stable  will  not 
turn  as  tightly  as  a  somewhat  less  stable  ship  so  that  a 
highly  stable  ship  may  compromise  maneuverability.   There- 
fore, a  ship  should  be  designed  for  a  moderate  amount  of 
stability  so  as  to  be  able  to  go  straight  without  compro- 
mising maneuverability  too  much. 

With  these  factors  in  mind  then,  the  stability  of  five 
different  ships  was  studied  in  Part  III  and  maneuverability 
in  Part  IV.   The  first  objective  of  this  thesis  was  to 
investigate  the  stability  of  a  ship  with  movable  fin.   The 
second  objective  was  to  investigate  the  maneuverability  of 
an  unstable  ship  in  standard  maneuvers  to  compare  with  those 
of  a  stable  ship.   A  third  objective  was  to  apply  DSL/360 
Language  to  simulate  the  ship  motions. 
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II.   EQUATION  OF  SHIP  MOTION  IN  THE  HORIZONTAL  PLANE 

As  developed  in  Ref.  2,  the  general  motion  of  a  ship  is 
that  of  a  rigid  body  in  six  degrees  of  freedom  subject  to 
gravity  and  buoyancy  forces  as  well  as  controlling  forces 
and  hydrodynamic  reactions  and  subject  to  hydrodynamic  or 
other  disturbances  or  excitations. 

Considering  axes  fixed  to  the  body  and  parallel  with 
the  principal  axes  of  inertia  but  with  arbitrary  origin  as 
shown  in  Fig.  1,  the  linear  and  angular  momentum  equations 
may  be  written  as 


where 


(1) 


m  =  mass  of  ship 

— *• 

U  =  linear  velocity  vector 

=  ML  +   v<|  +  w  k 

u  =  rate  of  surging 
v  =  rate  of  swaying 
r  =  rate  of  heaving 

=  angular  velocity  vector 

=  pT  +  q£  +  r£  (2) 

p  =  rate  of  roll  =  0 
q  =  rate  of  yaw  =  9 
r  =  rate  of  pitch  =  y 
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F  = 

X  = 
Y  = 


Z 

% 


K 

M 
N 
H 


1,1,1 

x    y    z 


force  vector  acting  on  the  body 

XT  +  Yj  +  Z~£  (3) 

hydrodynamic  force  on  body  along  x-axis 

hydrodynamic  force  on  body  along  y-axis 

hydrodynamic  force  on  body  along  z-axis 

moment  vector  acting  on  body  about  a 

fixed  point 

kT  +  uj+   N~£  (4) 

rolling  moment  about  x-axis 

pitching  moment  about  the  y-axis 

yawing  moment  about  the  z-axis 

angular  moment  vector 

I*p~i   +    tyqj+    IjrT  (5) 

mass   moment    of    inertia   about    the    axis 

x,    y   and   z   respectively. 


FIGURE    1.       Body    axes 
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Inserting  Equations  (2),  (3),  (4)  and  (5)  in  Eq .  (1) 
gives : 

y  -7n\v+r<*--pi>r-$*Cr\l?)+z*(fr-p)+Zt(tp*r)}      \(e) 
m  =z^+(ix-lt)rp+}n{^(^far--tr)-^(ur+pv--f^}   \(7) 

In  Eq .  (6)  and  (7),  the  left  hand  side  represents  the 
forces  and  moments  acting  on  the  ship  and  the  right  hand 
side  represents  the  dynamic  responses  of  the  ship, 
respectively.   The  forces  and  moments  acting  on  the  ship 
depend  on  the  geometry  of  the  hull,  the  motion  and  the  fluid, 
That  is, 

F  (or  /7l   )  =  f  (properties  of  ship,  properties  of 
fluid,  properties  of  ship  motion) 

Since  the  properties  of  ship  and  fluid  are  constant  for 
a  given  ship  in  a  given  fluid  without  excitation  forces,  the 
forces  or  moments  may  be  expressed  as  a  function  of  the 
motion  as, 

F(VL  Ik)  =  f  ("*>,  %  ^  ft  B/  ff  UL^WjP^jr 

(  8) 

°-,  *,  ^>;  A  h  r,  S/  S,  I  Jc) 
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where  x  ,  y   and  z   represent  a  system  of  reference  axes 
o    o      o 

fixed  in  the  earth  and 


/Y  H      ~sr     &    9    y     ~   orientation  parameters 

"/  v,  m  r,  h  r  i 

0-,  v-,  tir,  r,%,r\=  motion  para" eters 

Sji/j       -  contrd'  s„rf  ace  palters 

Considering  motion  in  unrestricted  calm  water,  only  three 
degrees  of  freedom  are  of  concern:  yaw,  sway  and  surge. 
Therefore,  the  forces  and  moments  are  independent  of  the 
change  in  orientation  with  respect  to  the  earth  axes. 

Furthermore  since  the  forces  and  moments  produced  on  the 
ship  as  a  result  of  q   and  £  are  normally  negligible. 

F(or/tt)  =  f(u,  v,  r,  u,  v,  r,  S)  (9) 

Since  the  motion  is  restrained  to  the  horizontal  plane 
for  the  surface  ship  in  calm  water  as  prescribed  in  the 
previous  section,  roll,  pitching  and  heaving  are  taken  as 
zero,  i.e.,  p=q=w=0  in  Eqs .  (6)  and  (7).   Furthermore,  most 
ships  have  q^~   0  so  that  Eqs.  (6)  and  (7)  may  be  written  as, 

x  =  TnftJL  -rir~r^rx) 

Y  =  #7(^+rU+/fcrf)  (10) 
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From  Eqs.  (9)  and  (10)  then,  these  become: 

x:  7r)(u-rv-x*f)  =  ktUyVsTj  u,v,r,&) 

Y:  m(v-^iiLY^r)  =.  \^Jv,r,u./vJrJ  6)       (11) 
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III.   DYNAMIC  STABILITY  OF  THE  SHIP 

Dynamic  stability  is  directly  related  to  the  magnitude 
of   yaw  and  sway  deviations  caused  by  small  initial  dis- 
turbances.  In  many  problems  related  to  dynamic  stability 
with  small  disturbances,  the  linear  analysis  of  motion  is 
quite  useful  and  will  be  applied  in  the  following  to  the 
present  problem. 

A.   LINEAR  EQUATION  OF  SHIP  MOTION 

To  simplify  the  equation  of  ship  motion,  the  following 
assumptions  are  made: 

(1)  The  ship  has  adequate  stability. 

(2)  The  water  is  so  calm  that  the  rolling,  pitching 
and  heaving  motions  are  negligible.   Therefore,  the  motion 
of  ship  may  be  regared  as  taking  place  in  the  horizontal 
plane  only. 

(3)  Yaw  and  sway  do  not  affect  the  forward  speed 
appreciably,  and  small  changes  in  forward  speed  do  not  affect 
yaw  and  sway  motions  so  that  the  surge  equation  is  decoupled. 

(4)  The  ship  has  one  rudder  at  the  stern  along  the 
centerline  and  the  effect  of  the  propeller  is  disregarded. 

(5)  The  forces  and  moments  on  hull  and  rudder  are 
expressed  by  means  of  velocity  and  acceleration  derivatives. 

In  order  to  linearize  the  equations  of  ship  motion  about 
an  initial  equilibrium • condition  of  motion  with  the  previous 
assumptions,  only  the  linear  terms  in  the  change  of  the  value 
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of  the  variables  from  equilibrium,  i.e.  Au,  Av,  Ar,  Au,  av, 
Ar  and  a&  are  considered. 
From  Eq.  (9)  then, 


X 


=  Xo+'XuaU+XvaW  Xr*r+  Xa*&  +X*av+Xs*£ 


Since  in  initial  equilibrium  V*  =  To  "=  uc  ~  VI  —  To  —  do  ~  0  . 

ait=  V-Vo  =  v' 

Ar -  r-To  =  r 

au  =■  u  -  do  --•  u 

Ar  -r  -h  "T 

Therefore , 
and  likewise, 


3) 


where  X  ,  Y   and  N   are  the  values  of  the  initial  equilibrium 
o '   o      o 

conditions. 

In  the  initial  equiliburium  conditions  of  straight  ahead 
motion  at  constant  speed,  there  are  no  forces  acting  on  the 
ship.   Hence 

X.  =  Y.  -  a/.  =  o 
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Furthermore,  since  the  ship  is  symmetric  about  the  x-z  plane 

X  =  /"  =•  Xu-=  Xr=  Xf  =  A/«  -  Na  =  o 

and 
Therefore, 

y  =  y^+yrr  +  /^i/fr-\-/r^  d4) 

A/  =  A/U+  A/r  I"  +  Viz  V  +  A/r  r  +  ^  <? 

Linearizing  Eq .  (10)  gives 

This  is  simplified  further  by  neglecting  second  order  terms 
so  that 


and  likewise, 

/^/=  T^  f  +  ^%  (v  +  ru0) 

From  Eqs.  (14)  and  (15) 
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(15) 


(16) 


For  convenience,  Eq .  (16)  is  non-dimensionalized .   Every 
term  of  the  first  and  second  equations  for  the  X-force  and 
the  Y-force  in  Eq .  (16)  has  the  dimensions  of  force  and  the 
third  equation  for  moments  has  the  dimension  of  a  moment. 

To  non-dimensionalize  these,  the  force  equations  are 

f     x 
divided  through  by  (-^rLnU)  and  the  moment  equation  is  divided 

9     2       2 
by  (~zl-nU)   where  f  =   dinsity  of  ship,  L  =  ship's  length, 

U  =  ship  speed  and  H  =  draft  of  ship.   This  gives, 


where 


T  -  J   /? 


for  small  perturbations, 


'    w  A 
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In  Eq.  (17)  the  dimensionless  time  is  defined  as 

U 


f=  +  t 


U 


B.   DYNAMIC  STABILITY  WITH  FIXED  CONTROL 

Taking  the  Laplace  Transform  of  Eq.  (17)  with  the  rudder 
fixed  at  6=  0  gi1 


gives : 


[M*  -m  tys+ a/v']v'(s)  +[MrT/)s  +  Mr!*  4  V)J  u*>o 


(18) 


Then,  using  Cramer's  rule 


v(s) 


r(s)  = 


0 


A 

0 


A 


where 


A 


Expanding  this  determinant  gives  the  characteristic 
equation 


A  ^A  s2+  B5  +  C 


(19) 


22 


where 


A  m(yt-m')(dt-&)-M-/ite)(fiJir»i%) 


rfr'(y/-»  '**)  -  (  Yr  m ')(Mr'-*>  %  )  (20 ) 


c  =  //(t/ra'Xf'J-A/vty-to') 


The    roots   of   Eq.     (19)    are 


(21) 


U  2      —  ~ 

2 

Therefore,  the  free  responses  of  ship,  i.e.,  sway  and 
yaw  motions  are  given  by 


<rrt  <r2t 


/ 


v(f)  =  v,e    +  vxe  (22) 

r  (f;  =^e     +  rze 


where  v.,  v~  ,  r..  and  r„  are  arbitrary  constants  depending  on 
initial  conditions. 

As  shown  in  Eq.  (22),  the  dynamical  stability  depends  on 
the  roots  of  the  characteristic  equation.   Since  (J]~ )>  (Fz 

Msfn.    v'(i')-0       OAid     <lwn  r'(t')=0 
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if  and  only  if    (f\   <  0 

If  U^yo   then,  y~(oo)=£  0      and  ^(00)^0     which  means  that  the  ship 
is  dynamically  unstable  in  straight  line  motion.   The  root, 
6]   is  called  the  'stability  index'. 

In  Appendix  B  and  computer  program  1,  a  five  foot  Series 
60  model  of  the  Davidson  Laboratory  with  block  coefficient, 
CL.  =  0.7,  and  no  propeller,  was  modified  to  be  (a)  stable 
(b)  marginally  unstable  (c)  unstable  (d)  very  unstable  and 
(e)  extremely  unstable.   These  models  are  referred  to  as 
model  A,  B,  C,  D  and  E,  respectively. 

The  hydrodynamic  derivatives  and  stability  indices  of 
these  models  are  in  Table  1. 

Eqs.  (21)  and  (22)  were  simulated  to  obtain  the  response 
of  each  model  to  small  a  disturbance.  (See  computer  program 
7  and  Fig.  2) 

As  shown  in  Fig.  2  (which  is  the  result  of  computer 
program  7),  the  yaw  rate,  r,  as  well  as  the  sway  velocity, 
v,  of  the  stable  model  A  decreases  as  time  increases  but 
the  unstable  models  B,  C,  D  and  E  tend  to  diverge  due  to 
the  positive  stability  index. 


TABLE  I.   HYDRODYNAMIC  DERIVATIVES  AND 
STABILITY  INDICES  OF  THE  FIVE  MODELS. 


mod- 
el 

/; 

ti 

// 

a// 

y> 

V* 

Yr 

a// 

07 

A 

-.3T44 

-.ofo2 

+.0774 

-,^3r 

-.mo 

+  .004-0 

-+-  .ooq.o 

-.0/32 

-.321*3 

B 

-•3/3/ 

-•//45" 

f.0^7 

-.  o$32 

->!U2 

+.0034. 

+.00^4 

-M*!\ 

+  .42^? 

C 

-.2/3/ 

-.fzU 

+  .044-6 

-.^4-35" 

-./W8 

+  .ooz4- 

-\-.ooz4- 

-.012-1 

+  .2^1 

D 

-»-J/3' 

-■144-1 

+.026$ 

-o3*d 

-.126% 

4-.  00/0 

+  .  0o\D 

-.0//7 

+.^94- 

E 

-.i%6j 

-■1717 

-  .0bb$ 

-.021 i 

-.iftoo 

0. 

0. 

-.ollZ 

4/.  34-3? 
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C.   STABILIZATION  OF  THE  DYNAMICALLY  UNSTABLE  SHIP  WITH 

FIXED  RUDDER 

As  described  in  the  previous  section,  a  dynamically 
unstable  ship  has 


c   -f+/y-^>   >0 


To  stabilize  this  unstable  ship,  (ff  must  be  made  less 
than  or  equal  to  zero,  i.e.,  d?^0  or, 

B 


fiP+W)  <  -* 


By  the  nature  of  hydrodynamic  derivatives  for  the  ship 
A^  0  and  B\  0  always.   Therefore,  C  must  be  greater  than  or 
equal  to  zero  for  stability.   Thus,  C  is  called  the  "stability 
criterion . " 

From  Eq.  (20) 

C  =X/(a//-  m'^')-^(y;-  m')^0         (20a) 


Then,  considering  the  characteristics  of  each  of  the 
hydrodynamic  derivative,  mass  and  the  distance  of  the  C.G 
from  the  origin  of  the  fixed  coordinates  on  the  ship. 

Y  is  always  large  negative  quantity 
N  is  always  large  negative  quantity 
N  may  be  positive  or  negative 

Y  is  small  but  may  be  positive  or  negative 
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m    is  always  large  and  positive 

X„   is  always  small  but  may  be  positive  or 

negative 

Finally,  N   has  the  greatest  effect  on  dynamical  stability 

since  the  term  (/r~^)  is  a  large  negative  quantity  in 

equation  (20a).   If,  for  example,  N   becomes  more  positive, 

the  stability  criterion  C  will  become  more  positive  so  that 

the  stability  of  the  ship  will  be  improved.   It  is  noted 

that  N   can  be  controlled  by  a  lifting  surface. 

1 .   Sensitivity  of  Dynamic  Stability  to  a  Lifting 

Surface 

Equation  (21)  was  simulated  using  the  computer 

program  DSL/360  Language  for  the  effect  of  the  location,  size 

and  aspect  ratio  of  a  lifting  surface  (movable  fin)  on  the 

models  A,  B,  C,  D  and  E.  (See  the  computer  program  2).   The 

results  are  summarized  in  Table  II  and  III. 

The  following  are  indicated  by  Tables  II  and  III  and  by 

Figs.  3,  4  and  5: 

1)  The  position  of  the  fin  at  which  the  stability  index 

is  maximum  is  independent  of  the  area  and  aspect  ratio  of  the 

fin,  because  this  position  is  the  pivot  point  at  which  the 

relative  flow  is  parallel  to  the  centerline  of  the  ship. 

(See  Fig.  14) 

The  parallel  flow  makes  no  attack  angle  on  the  fin  so 

that  Y  ,  N  ,  Y   and  N   of  the  fin  are  zero  but  Y*   N' ,  Y*  and 
vvr      r  vvr 

N-  are  not  zero  as  described  in  Appendix  A.   Therefore,  the 
stability  criterion  C  in  Eq.  (20)  is  the  same  as  that  of  the 
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TABLE  II.  EFFECT  OF  FIN  AREA  AND  FIN  LOCATION  WITH 
FIXED  ASPECT  RATIO  ON  DYNAMICAL  STABILITY 


Fin 

Xf/L 

A 

B 

C 

D 

E 

No  fin 

-0.3213 

+0.0238 

+0.2561 

+0.5194 

+1.3438 

OL   ..=1.5 

eff 
Af =0.0056 

hf=0.015 

-  0.5 

-  .4366 

-  .1066 

+  .1108 

+  .3625 

+1.1684 

m 

-  .3492 

-  .0020 

-  .2324 

+  .4968 

+1.3223 

pivot 
point 

-  .32079 

+  .0242 

'+  .2569 

+  .5197 

+1.3441 

(0.48) 

(0.4) 

(0.34) 

(1.305) 

(0.27) 

+  0.5 

-  .3208 

+  .0226 

+  .2512 

+  .5109 

+1.3288 

#   =1.5 
eff 

A  =. 01265 
f 

hf=.0225 

-  0.5 

-  .5707 

-  .2616 

-  .0642 

+  .1710 

+  .9474 

& 

-  .3834 

-  .0341 

+  .2026 

+  .4682 

+1.2944 

pivot 
point 

-  .3203 

+  .0246 

+  .2567 

+  .5199 

+1.3442 

(0.48) 

(0.4) 

(0.34) 

(0.305) 

(0.27) 

+  0.5 

-  .3204 

+  .0212 

+  .2458 

+  .5011 

+1.3110 

CL   ..=1.5 
eff 

Af=.0225 
hf=.03 

-  0.5 

-  .7376 

-  .4583 

-  .2887 

-  .0780 

-  .6717 

£ 

-  .4299 

-  .0780 

+  .1617 

+  .4286 

+1.2548 

pivot 
point 

-  .3197 

+  .0250 

+  .2571 

+  .5202 

+1.3444 

. (0.48) 

(0.4) 

(0.34) 

(0.305) 

(0.27) 

+  0.5 

-  .3199 

+  .0196 

+  .2393 

+  .4893 

+1.2888 

d   ..=1.5 
eff 

Af =0.0506 
hf=0.045 

-  0.5 

-1.0956 

-  .8930 

-  .7921 

-  .6486 

-  .0573 

& 

-  .5532 

-  .1963 

+  .0507 

-  .3198 

+1.1427 

pivot 
point 

-  .3185 

+  .0260 

+  .2579 

+  .5209 

+1.3448 

(0.48) 

(0.4) 

(0.34) 

(0.305) 

(0.27) 

+  0.5 

-  .3188 

+  .0162 

+  .2259 

+  .4643 

+1.2393 

destabiliz 
area 

ing 

0.415i-^ 
0.5L 

0.35L~ 
0.445L 

0.3L  ^- 
0.375L 

0.27L  — 
0.34L 

0.24L~ 
0.3L 

Note:  All  values 
in  (  )  wh 
(T*  of  mode 


in  column  A,B,C,D  and  E  (except  the  values 
ich  mean  pivot  point  distance  from  20  are 
1  A,B,C,D  and  E,  respectively. 
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TABLE  III.   EFFECT  OF  FIN  LOCATION  AND  ASPECT  RATIO  WITH 
FIXED  AREA  ON  DYNAMICAL  STABILITY 


Fin  Type 

Xf/L 

Model  A 

B 

C 

D 

E 

No  fin 

-0.3213 

+0.0238 

+0.2561 

+0.5194 

+1.3438 

%=    1 .  5 
Af =0.0056 

hf=0.015 

-  0.5 

-0.4366 

-0.1066 

+0. 1108 

+0.3625 

+1. 1684 

X 

-0.3492 

-0.0020 

+0.2324 

+0.4968 

+1.3223 

pivot 
point 

-0.3208 

+0.0242 

+0.2569 

+0.5197 

+1.3441 

(0.48) 

(0.40) 

(0.34) 

(0.305) 

(0.27) 

+  0.5 

-0. 3208 

+0.0226 

+0.2512 

+0.5109 

+1.3288 

*%=  2. 

Af =0.0056 

hf=0.0173 

-  0.5 

-0.4557 

-0. 1275 

+0.0872 

+0.3373 

+1.1442 

& 

-0.3539 

-0.0063 

+0.2285 

+0.4932 

+1.3189 

pivot 

-0.3207 

+0.0242 

+0.2565 

+0.5197 

+1.3439 

point 

(0.48) 

(0.4) 

(0.34) 

(0.305) 

(0.27) 

+  0.5 

+0.3208 

+0.0224 

+0.2505 

+0.0510 

+1.3265 

Af=0.0056 
hf=0.0212 

-  0.5 

-0.4823 

-0.1572 

+0.0543 

+0.3021 

+1.1036 

® 

-0.3605 

-0.0122 

+0.2231 

+0.4881 

+1.3164 

pivot 
point 

-0.3206 

+0.0242 

+0.2565 

+0.5197 

+1.3440 

(0.48) 

(0.40) 

(0.34) 

(0.305) 

(0.27) 

+  0.5 

-0.3207 

+0.0221 

+0.2495 

+0.5078 

+1.3235 

%=4. 

Af =0.0056 

hf=0. 02448 

-  0.5 

-0.5000 

-0.1769 

+0.0325 

+0.2788 

+1.0786 

$ 

-0.3649 

-0.0162 

+0.2195 

+0.4847 

+1.3112 

pivot 
point 

-0.3206 

+0.0244 

+0.2565 

+0.5198 

+1.3441 

(0.48) 

(0.4) 

(0.34) 

(0.305) 

(0.27) 

+  0.5 

-0.3206 

+0.0219 

+0.2488 

+0.5066 

+1.3215 

Note:  All  values  in  column  A,B,C,D 
in  (  )  which  mean  pivot  poi 
(J*r  of  model  A,B,C,D  and  E,  r 
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and  E  (except  the  values 
nt  distance  from  35.)  are 
espectively . 


iroiel  E 


1   rrodel  A. 


.022£ 
.0505 


.022£ 

.030 

.01£ 


-.5 


Fig 


-.3 


.1  0.1  .3 

fin  position 


3   The  effect  of  fin  location  and  fin  area  with 
constant  aspect  ratio  on  stability  of  the  ship 
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model  E 


75 


curve  effective 

from  top     aspect  ratio 


1st 

l.S 

2nd 

2.0 

3rd 

3.0 

Uth 
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_4   model  D 


a   model  G 


_2    model   B 


model  A 


*fin  area;      *  .00£5 

L _  ; 

1         .3         .5 
fin  position  Xf  /  L 

FIGURE  4.   The  effect  of  fin  location  and  aspect  ratio  with 
constant  area  of  fin  on  the  stability  of  the  ship. 
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no-fin  condition.   In  other  words,  the  stability  criterion, 
C  is  minimum  at  the  pivot  point  during  shifting  of  the  fin 
from  stern  to  bow. 

This  pivot  point  is  the  position  of  the  fin  at  which  it 
exerts  maximum  destabilizing  influences  on  the  ship. 

2)  The  effect  of  the  fin  on  stability  increases  with 
fin  size. 

3)  The  more  unstable  the  ship,  the  closer  the  pivot 
point  moves  toward  the  midship. 

4)  Compared  to  the  no-fin  condition,  the  placement  of 
the  fin  near  the  pivot  point  has  a  negative  effect  on 
stability.   This  criterion  is  called  the  "destabilizing  area" 
in  this  thesis.   This  destabilizing  area  is  due  to  the 
relative  change  in  A,  B,  and  C  in  Eq.  (21)  as  the  fin  location 
is  shifted.   C  is  almost  constant  near  the  pivot  point  and 
the  destabilizing  area  is  independent  of  the  size  and  type 

of  the  fin. 

5)  Moving  the  fin  toward  the  stern  from  midship 
improves  the  dynamic  stability  of  the  unstable  ship.   Con- 
versely, moving  the  fin  forward  of  midship  does  not  cause 
much  change  from  no-fin  condition  regardless  of  the  size  or 
type  of  the  fin. 

6)  Stability  is  sensitive  to  the  fin  area  but  rather 
insensitive  to  changes  in  the  aspect  ratio  of  the  fin. 

2 .   Sensitivity  of  Stability  of  the  Unstable  Ship  to  the 
Position  of  the  C.G.  of  the  Ship 

It  is  clear  that  moving  the  C.G.  toward  the  bow  tends 
to  make  'C  more  positive  in  Eq .  (20).   The  effect  of  moving 
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C.G.  was  studied  by  use  of  Eqs .  (20)  and  (21)  using  DSL/360 
Language  for  model  C.   The  results  are  shown  in  Table  IV. 

(See  computer  program  3) . 

i 

TABLE  IV.   THE  EFFECT  OF  LOCATION  OF  C.G.  ON 
THE  STABILITY  INDEX 


XG/L 

-0.4 

-0.25 

X 

+0.05 

+0.1 

+  0.25 

<Ti 

+    .6424 

+    .4022 

+    .0775 

+    .015 

-    .0473 

-    .2297 

As  shown  in  the  above  table,  stability  was  improved  by 
moving  the  C.G.  toward  the  bow  as  expected.  However,  this 
is  not  a  practical  method  of  stabilizing  an  unstable  ship. 


D.   STABILIZATION  OF  THE  DYNAMICALLY  UNSTABLE  SHIP  WITH 

RUDDER  CONTROL 

In  previous  sections,  the  stabilization  of  a  dynamically 
unstable  ship  using  lifting  surface  was  discussed.   However, 
it  is  clear  that  the  dynamically  unstable  ship  can  be  made 
directionally  stable  if  the  rudder  can  be  forced  to  compen- 
sate for  disturbances  so  that  the  ship  will  maintain  course 
with  only  small  oscillations  in  yaw  and  sway.   If  the 
deflection  of  the  rudder,  £  ,  is  proportional  to  the  heading 
error,  y",  and  the  angular  velocity,  r  =  ^1/ ,    then  the  equation 
of  rudder  deflection  becomes 


&({)  =  kia'i  +  kisd') 


(23) 
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where/?.  and/2„  are  constants  of  proportionality  of  the 
control  system. 

Inserting  Eq .  (23)  into  Eq .  (17)  and  taking  Laplace 
Transform  gives 

(24) 

[(ti}-to%)s+fil^v(s)+lfc/t  -l/)S  V  Mt'-&  %'+k  d/js+k  tJi)f(s)  =  o 


As  these  algebraic  equations  of  (24)  are  homogeneous,  a 
non-trivial  solution  (v  %    y   ^  0)  demands  the  determinant  of 
the  system,  i.e.,  characteristic  equation,  to  be  zero. 
Let  the  characteristic  equation  be 

As'+Bs'+Cs  +  i)  =  o  (25) 


where 


a/;  (# '-  w  'V;  -  c  Yt'-M  '■+  £  #W~*  '^' } 
a/;  ( y,  -  » + 4  // )  -  *.  #'m'  -*»  V  J 

p  =  k,  y^i  -  k  v;  yI 
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/    / 


Let  the  roots  of  Eq.  (25)  be   (Jf  (f^   and  Qz' .      Then, 


j 


fu'i^e^te^^i1' 


«-      Xf     .jr,v  <26> 


T/?ifV=  i^e  +ijg    + 1^<? 


^ 


If  and  only  if  CTT ,  CT\.  and  (J3  are  negative,  then  y(-oo)  and 
V(oo)  will  be  zero  and  the  unstable  ship  will  be  stable. 

Furthermore,  if  the  roots  are  three  repeated  roots,  i.e., 

*     /     / 

(Tf  =  (Ti  =  (Tj  <^  0,  the  ship  will  be  stable  with  the  shortest 

possible  response. 

By  designing  for  the  optimum  value  of  k,  and  kz  in  the 
control  mechanism,  the  ship  will  be  automatically  controlled 
for  suitable  stability. 

To  find  the  three  repeated  negative  roots  of  the 
characteristic  Eq.  (25),  the  'MITROVIC  METHOD'  was  applied. 
(See  Appendix  C). 

In  computer  program  4,  the  optimum  K,   and  r?z  for  the  three 
repeated  negative  roots  of  Eq .  (25)  were  calculated  for  model 
C  and  Eq .  (26)  was  simulated  by  the  DSL/360  Language.  (See 
computer  program  4). 

As  shown  in  Fig.  5-1  which  is  the  result  of  computer 
program  4,  the  response  with  proper  /v,  and  az  for  three 
negative  real  roots  of  characteristic  equation  (25)  is  an 
exponential  transient.   Conversely,  if  the  auto-steering 
device  has  some  kt   and  fiz  which  give  complex  roots  of  charac- 
teristic equation,  the  response  will  be  stable  with  some 

oscillations . 
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The  auto-steering  device  should  be  able  to  select  the 
proper  control  costants,  ki   and  /?x  . 

It  is  evident  that  the  dynamically  unstable  ship  can 
be  made  to  be  stable  directionally  with  auto-steering  control 
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IV.   MANEUVERABILITY  OF  THE  UNSTABLE  SHIP 

To  predict  the  performance  of  a  dynamically  unstable  ship 
as  well  as  to  get  realistic  predictions  of  maneuvers  with 
large  rudder  angles,  a  non  linear  mathematical  model  is 
needed. 

A.   NON-LINEAR  EQUATION  OF  SHIP  MOTION. 

In  the  non-linear  analysis  of  the  ship  motion,  the 
following  assumptions  which  are  rather  typical  in  the 
literatures  Refs.  2  and  8  were  applied. 

1)  Terms  of  higher  order  than  third  order  are  neglected. 

2)  Rudder  force  and  moment  derivatives  of  higher  than 
first  order,  and  effects  of  rudder  angular  rate  are  negligible 

3)  Odd  terms  in  v,  r,  <£  ,  v  and  r  in  the  X-force  equation 
and  even  terms  in  the  same  parameters  of  the  Y-force  and 

N-moment  equations  are  eliminated  due  to  symmetry. 

4)  Cross  coupling  terms  between  acceleration  and 
velocity  are  neglected.   Taking  the  Taylor  Series  expansion 
of  Eq .  (11)  with  the  above  assumptions: 

U   z:   - , 
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r  '=  (l»'-YM(«Wy  h  S)  -  (■» %'-*/*) i «  v,  h  <0 


where 


f  ?> 


fyWhl)-  A//+A/0«4u'+/iuu4u'  +  A/„'v^O'4 


_1      ;/ 


(28) 
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Integrating   Eq .    (27) 

j.  * 

u'(i')=f  u'(t')  di' 

(*' 

v'(i')  -  J   v-'(f)dt 

ft'  (29) 

o 

Note:   The  time  domain  should  be  non-dimensional. 

So  far  the  reference  axes  x,  y  and  z  are  fixed  to  the 

moving  ship.   But  it  is  convenient  to  take  the  reference  axes 

x   and  y   fixed  to  the  earth  in  order  to  determine  the  path 
o      J  o  ^ 

of  the  ship.   Adopting  the  figure  from  Ref.  1  gives 


->  2-c 


V 


FIGURE  6.   Relations  between  axes  fixed 
in  the  earth  and  axes  fixed  in  the  ship. 
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The  following  relationships  may  be  written  as 

%ii)~\x  (t)  cm  i>(t)  -  v(t)  /tin  f  (t) 

(30) 

yo  (t)  =  a  (t)/iu\  fit)  -f-  vlt)  cm.  fit) 


'a 


(31) 


Non-dimensionalizing  Eqs .  (30)  and  (31)  using  the 
definitions. 


*'-r  ,•  *'-£>£/- 

and, 

Z-fe*'&T* 

-     u 


or 


X'tf',=  \*'A* 


(32) 


Likewise , 


i'tf)=)1  Jt 
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B.   TURNING  CIRCLE  WITH  SELECTED  RUDDER  DEFLECTION 

Fig.  7  shows  the  typical  rudder  deflection  with  respect 
to  time  in  a  turning  maneuver. 


& 


MM 


S  *w 


Z^Tte 


FIGURE  7.   Rudder  deflection  vs 
time  in  turning  maneuver 


Using  the  ramp  function  for  rudder  deflection,  then 


J (i)  =  o(  [  rampd-fa*)  -  rmplt-t^} 


(33) 


where 

O   =  rudder  deflection 
ramp  (t)  =  ramp  function 

C\  =  rudder  deflection  rate 
*j£ojj  =  time  delay 

<^<W=  time  at   c(>r>a* 
Calculation  of  steady  turning  radius. 

The  relation  of  the  angular  velocity  with  linear 

velocity  gives: 

(radius)  x  (angular  velocity)  =  (linear  velocity) 
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or     Radius 


U(t) 

rep 


U.\t)  +  V(t) 


ret) 


Taking  non-dimensional  terms 


L 


u     ilu^v7- 


U  L 


_  J(u/uT+(Vu): 


r 


r(-£) 


or 


RaJ   - 


J(u?+(V) 


(34) 


where 


Rod     = 


K  a  Jili 


fUL 


L- 


Figure  8  shows  the  typical  turning  trajectory  and 
terminologies  corresponding  to  Table  5.   Eqs .  (31)  and  (34) 
were  simulated  for  the  5  models  using  computer  program  5  in 
turning  circle  maneuvers  with  a  rudder  deflection  rate  of 
2  degrees  per  second  for  the  prototype  (L  =  500')  with  zero 
time  lag  and  constant  final  angles  of  10  ,  20   and  35  . 
The  results  are  presented  in  Table  5  and  Figs.  8a,  8b  and  8c 
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Definition  of  turning  circle  parameters 
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FIGURE  8a.   Turning  circle  of  five  models  with 
35  degrees  of  rudder  angle. 
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FIGURE  8b.   Turning  circle  of  five  models  with 

25  degrees  of  rudder  angle. 
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FIGURE  8c.   Turning  circle  of  five  models  with 
10  degrees  of  rudder  angle. 
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TABLE  V.   CHARACTERISTICS  OF  MODEL  A,B,C,D  AND  E 

IN  TURNING  CIRCLE 


J 

s 
o 
a 
<x> 
(-■ 

Adv. 

Max. 
Adv. 

Tran . 

T.D. 

Max . 
Tran. 

S.T.R. 

'  rev , 

10° 

A 

6.86 

6. 8825 

4.86 

10.30 

10.322 

5.0722 

31.6 

B 

5.30 

5.3224 

3.25 

7.08 

7. 1079 

3.4498 

21.4 

C 

4.49 

4.5088 

2.65 

5.88 

5.7079 

2.7611 

17.2 

D 

3.94 

3.9667 

2.18 

4.79 

4.8122 

2.3290 

14.6 

E 

3.08 

3.1277 

1.45 

3.51 

3.5490 

1.7295 

10.  8 

20° 

A 

4.80 

4 . 8269 

3.16 

6.91 

6.9411 

3.4177 

21.2 

B 

4.07 

4.0996 

2.40 

5.46 

5.4962 

2.6946 

16.8 

C 

3.60 

3.6306 

2.06 

4.66 

4.6863 

2.2951 

14.2 

D 

3.29 

3.3134 

1.80 

4.10 

4.1379 

2.0266 

12.6 

E 

2.70 

2.7892 

1.35 

3.23 

3.2719 

1.6070 

10.0 

35° 

A 

3.840 

3.8685 

2.33 

5.2203 

5.2433 

2.6364 

16.2 

B 

3.430 

3.4584 

1.98 

4.435 

4.4551 

2.1938 

13.6 

C 

3.  120 

3.1538 

1.77 

3.900 

3.9380 

1.9385 

12.0 

D 

2.918 

2.9475 

1.49 

3.540 

3.5804 

1.760 

10.8 

E 

2.550 

2.5943 

1.15 

2.790 

2.9890 

1.4727 

8.6 

Note:   Adv.  =  Advance/L 

Tran.  =  Transfer/L 

T.D.  =  Tactical  Diameter/L 

S.T.R.  =  Steady  turn  radius/L 


48 


C 

o 

•H 

TJ 

-P 

U 

f  } 

ri 

a> 

H 

w 

«m 

^ 

<D 

Q)  t3 

-P 

o 

U 

fc; 

0) 

<f 

T) 

U 

X) 

rt 

0 

a, 

^ 

0) 

<M 

H 

O 

f  > 

S< 

•H 

Q) 

o 

a; 

M 

fo 

£, 

<P 

•H 

xi 

>-; 

^  IT 

0 

o^ 

^ 

-P 

t 

ri 

G 

(!) 

X 

^: 

u 

-P 

Ti 

0) 

U 

-O 

•H 

to 

S>> 

C 

-P 

O 

•H 

■H 

rH 

4-> 

•H 

r< 

-D 

rH 

crt 

o 

-P 

r; 

W 

On 

• 

W> 

•ri 

[r, 

49 


w 
u 


bO 
U 


O 

u 

•H 
O 

bO 

c 

•H 

c 


•H 

3 


to 
o 

I 


bO 
•H 


VX> 

1/  sout-Apv 


50 


•n     <        cq 
g                / 

OQ            K 

//I] 

/. 

///]// 

MM 

7// 

1 
1 

/ 

/      / 

(A 

a) 

H 

o 

•H 
O 

- 

•H 

(D 

fc 

i-t 

^ 

W) 

3 

£ 

-P 

ej 

C 

^ 

•r-1 

C) 

TJ 

<D 

TJ 

(H 

2 

bO 

k 

<L> 

IS 
2 

ft! 

'M 

• 
to 

> 

u 

0) 

-p 

o 
•H 

+^ 

o 

o 

rH 

H 

•H 

51 


w 

U 

0) 


o 

to 


C 


W 

to 

•H 
T5 

to 
c 

•n 
C 

0) 

-p 

en 

CM 


•H 


00 


52 


w 

a) 

u 


W) 

c 

u 

Q> 


rH 
O 

•H 
O 

M 

•H 

C 

•H 


t>0 


5 

I 

•H 
-P 

tx 
c 

•H 

H 

B 

•H 

ri 

f\ 

rH 

• 

•H 
tt, 


(sseioT9J8iixp)^xi  OMxICHiC 


53 


Discussion  of  turning  circle  maneuver 

1)  As  shown  in  Fig.  9,  the  dynamically  less  stable  ship 
has  a  smaller  tactical  diameter,  advance,  steady  turning 
radius  and  circling  time  than  the  more  stable  ship,  because 
the  yawing  rate  of  the  unstable  ship  is  greater  than  that  of 
stable  ship  at  the  same  rudder  deflection.   However,  the 
turning  circle  parameters  are  becoming  constant  as  the 
stability  index  increases. 

2)  As  shown  in  Figs.  10,  11,  12  and  13,  the  rate  of 
change  of  the  turning  circle  parameters  versus  the  rudder 
deflection  is  high  at  the  small  rudder  angle.   In  other 
words,  the  above  parameters  are  less  sensitive  to  rudder 
angle  as  the  rudder  deflection  increases. 

3)  The  more  unstable  the  ship,  the  less  sensitive  the 
above  parameters  are  to  the  rudder  deflection. 

1 .   The  Effect  of  the  Location  of  the  Rudder  on  the 

Turning  Circle 

A  change  in  the  location  of  the  rudder  has  multiple 
effects  on  the  ship  motion.   One  is  it's  effect  on  dynamical 
stability,  which  is  the  same  as  the  effect  of  a  movable  fin 
described  in  section  (II-C-1)  and  the  other  is  it's  effect 
on  maneuverability. 

As  described  in  Appendix  A  and  Ref.  1, 

1)  Y  ,  Y-  and  Y^  are  independent  of  the  position  of 
the  rudder. 

2)  N  ,  Y  ,  N*  and-Y-  are  less  positive  if  the  rudder 

v    r    v      r 

moves  forward  from  the  stern. 
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3)  N   and  N*  are  less  negative  if  the  rudder  moves  from 
the  stern  to  the  midship  position  and  less  positive  from 
midship  to  the  bow. 

4)  Since  the  relative  flow  passing  the  ship  depends  on 
the  position  considered,  the  attack  angle  depends  on  the 
location  of  the  rudder.   Forward  of  the  pivot  point  the 
attack  angle  is  larger  than  rudder  angle,  at  the  pivot 
point  the  attack  angle  is  equal  to  rudder  angle  and  aft  of 
the  pivot  point  the  attack  angle  is  smaller  than  the  rudder 
angle. 

5)  N(  strongly  depends  on  the  location  of  the  rudder 
since  N^  =  Y^  -"Xr  where  xf  is  the  distance  of  the  rudder 
position  from  midship  and  is  negative  aft  of  midship.   Nr 
is  less  negative  if  the  rudder  moves  forward  to  the  bow, 
but  Nee  is  less  positive  if  o  is  negative.   This  less 
positive  effect  decreases  the  yawing  rate,  r,  due  to 
decreases  of  f„  in  Eqs.  (23)  and  (28).   Decreasing  r  in 
turn  causes  the  turning  radius  to  increases. 

On  the  other  hand,  the  movable  fin  effect  of  shifting 
the  rudder  forward  makes  the  ship  more  unstable  which  makes 
the  turning  radius  decrease.   These  two  compromising  effects 
are  shown  in  Fig.  14  and  15  which  were  obtained  from 
computer  program  6.   In  computer  program  6,  model  C  was 
simulated  for  the  evaluation  of  the  effect  of  the  rudder 
position  on  turning  maneuvers. 

As  shown  on  Fig.  14  which  represents  the  effect  of  the 
rudder  position  on  the  turning  radius  without  replacement  of 
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TABLE  VI.   THE  EFFECT  OF  RUDDER  POSITION  ON  THE  TURNING 
CIRCLE  OF  MODEL  C  WITHOUT  A  FIN  AT  THE  STERN. 


*f 

t  /AM/. 

1/ 

r 

Ra.Ji' 

€ 

-0.5 

7.  12 

-0.1617 

0.4413 

2.2960 

0.25616 

-0.4 

6.84 

-0.16571 

0.45704 

2.25 

-0.3 

6.72 

-0. 1701 

0.4688 

2.  164 

0.5482 

-0.22 

6.64 

-0.1735 

0.4734 

2.144 

0.6492 

-0.2 

6.6 

-0.1743 

0.4739 

2.14 

-0.17 

6.64 

-0.17556 

0.4742 

2.1412 

0.70581 

-0.10 

6.64 

-0.1784 

0.424 

2.15 

TL 

6.76 

-0. 1822 

0.4639 

2.19 

0.8552 

+  0.  1 

7.0 

-0.1856 

0.4485 

2.268 

0.90969 

0.15 

7.  16 

-0.1871 

0.4383 

2.32 

0.20 

7.44 

-0. 1885 

0.4265 

2.386 

0.9394 

0.21 

7.43 

-0.18869 

0.42398 

2.4 

0.9410 

0.22 

7.45 

-0.1889 

0.4214 

2.415 

0.9424 

0.24 

6.08 

+0.1345 

-0.51396 

1.963 

0.9445 

0.25 

6.08 

+0. 1352 

-0.5163 

1.955 

0.27 

6.04 

+0.1366 

-0.5205 

1.94 

0.30 

6.0 

+0.1385 

-0.5260 

1.92 

0.9456 

0.34 

5.92 

0. 14099 

-0.53204 

1.898 

0.9420 

0.5 

5.8 

0.1493 

-0.5425 

1.804 

0.8980 
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Mg.lU  The  effect  of  the  position  of  the  rudier  or 
turning  radius  without  a  fin  at  the  stern. 
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TABLE  Via.   THE  EFFECT  OF  RUDDER  POSITION  ON 
THE  TURNING  CIRCLE  WITH  A  STABILIZING  FIN  AT  THE  STERN 


"p '&<UfsCcru 

V 

r' 

KauL 

<K 

-0.5 

'       -0.1193 

0.3053 

3.299 

-0.3804 

-0.4 

-0. 1233 

0.3112 

3.238 

-0.2570 

-0.35 

-0. 1252 

0.3125 

3.225 

-0.1937 

-0.33 

-0.1259 

0.3127 

3.223 

-0.1683 

-0.3 

-0.1269 

0.3126 

3.225 

-0.1305 

-0.25 

-0.1286 

0.3113 

3.238 

-0.0625 

-0.2 

-0.  13 

0.3086 

3.268 

-0.0088 

-0.1 

-0.1322 

0.2983 

3.382 

+0.1001 

0. 

-0. 1333 

0.2810 

3.590 

+0.1892 

0.2 

-0. 1292 

0.2220 

4.542 

+0.2948 

0.24 

-0. 1267 

0.2054 

4.908 

+0.3039 

0.26 

-0.1251 

0. 1961 

5.  138 

+0.3070 

0.27 

-0. 1242 

0.1913 

5  .  268 

+0.308 

0.28 

+0.0572 

-0.2827 

3.544 

+0.309 

0.30 

+0.0621 

-0.2971 

3.373 

+0.310 

0.4 

+0.0779 

-0.3413 

2.939 

+0.305 

0.5 

+0.0874 

-0.3639 

2.758 

+0.2820 

-0. 1617 

0.4413 

2.296 

+0.2561 
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FIGURE  15.   The  effect  of  the  rudder  position  on  turning 
circle  of  model  C  with  a  stabilizing  fin  at  the  stern. 


59 


the  fin  at  the  original  rudder  location,  the  steady  turning 
radius  decreased  as  the  rudder  was  moved  from  -0 . 5L  to 
-0.2L,  and  increased  as  the  rudder  was  moved  from  -0 . 2L  to 
+0.22L.   At  the  rudder  position  0.24L,  the  turning  radius 
decreased  suddenly  from  2.415L  to  1.963L.   Also,  since  the 
turning  moment  at  this  point  was  reversed  in  direction,  the 
ship  suddenly  swings  to  the  opposite  side .   This  phenomenon 
is  the  characteristic  of  the  unstable  ship  as  described  in 
section  IV-D. 

If  the  turning  moment  is  decreased  due  to  shifting  the 
position  of  the  rudder,  the  yawing  rate  decreases  until  the 
rudder  change  occurs.   After  changing  the  turning  direction, 
the  magnitude  of  yaw  rate  was  suddenly  increased.   Finally, 
the  turning  radius  was  decreased  markedly. 

To  find  the  optimum  position  of  the  rudder  for  proper 
stability  and  maneuverability  of  model  C,  the  rudder  was 
moved  forward  to  the  bow  and  in  it's  location  a  stabilizing 
fin  was  installed. 

As  shown  in  Fig.  15  and  Table  Via,  the  turning  radius 
at  rudder  position  -0.33L  is  a  minimum  before  a  sudden 
change  of  rotation  and  in  addition  the  stability  index  is 
negative.   Although  this  minimum  turning  radius  (3.223L)  is 
greater  than  that  of  the  original  model  C  (2.296L),  it  is 
an  acceptable  optimum  rudder  position  due  to  it's  better 
stability  index,  i.e.,  negative  index. 
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C.   ZIG-ZAG  MANEUVER 

According  to  the  standard  zig-zag  maneuver,  the  relation 
of  the  rudder  deflection  (  d  )  and  ship's  heading  (y   )  is  as 


described  in  Fig.  16 
Sli) 

OwaX  ~ 


-i 


\r>cv,-\ 


FIGURE  16.   Relations  between  rudder  deflection, 
ship's  heading  and  time  in  standard  zig-zag  maneuver 


Using  a  ramp  function  for  the  rudder  displacement,  the 
rudder  angle  is  given  by: 

Sat)  =  o&roAup  (/-/; ;  -  mpit-L )  -  ra^p(t^H 

Tcmf  (x  -£t.)+  roAAf  ( p  Xs )  ] 


where 


C\    =  rudder  rate  (degree/sec.) 

Y0A(tf?(t)=  ramp  function 

/  =  time  delay  of  rudder 

/;  =  time  at  o  =  Cn-w 

/$  =  time  at  6  =  'f 

*4  -  time  at  6="crr-o-K 

ts=  time  at  I  =  "Jr 
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FIGURE  17.   Rudder  Deflection  and  Ship's  Heading  in  Z-maneuver 
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FIGURE  18.   Yawing  Rate  and  Sway  Velocity  in  Z-maneuver 
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curve  #1   :    path  of  model  A 
curve  #2    :   path  of  rodel  0 
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FIGURE  19.   Path  Trajectory  in  Z-maneuver. 
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In  computer  program  8,  the  standard  zig-zag  maneuvers 
were  simulated  for  model  A  and  C.   As  shown  in  Figs.  18,  19 
and  20  which  represent  the  results  of  the  program  8,  the 
following  can  be  observed: 

1)  The  amplitude  of  overshoot  in  yaw  angle  of  the 
unstable  ship  is  larger  than  that  of  the  stable  ship. 

2)  Although  the  first  reversing  execution  of  the  rudder 
of  the  unstable  ship  is  earlier  than  for  the  stable  ship, 
the  period  of  the  stable  ship  is  shorter  than  that  of  the 
unstable  ship. 

3)  Amplitude  of  overshoot  of  path  in  the  unstable  ship 
is  greater  than  stable  ship. 

4)  The  reason  for  1),  2)  and  3)  is  that  the  yaw  rate  of 
the  unstalbe  ship  is  greater  than  that  of  the  stable  ship  but 
the  sway  velocities  are  almost  the  same  for  either  ship  so 
that  it  takes  longer  for  the  unstable  ship  to  get  the  reverse 
heading  after  reversing  the  rudder  angle. 

D.   SPIRAL  MANEUVER 

The  standard  spiral  maneuver  of  model  C  was  carried  out 
using  DSL/360  in  computer  program  9.   The  process  was 
started  with  10  degrees  port  rudder  and  1  degree  of  deflect- 
ing increment  and  ended  with  10  degrees  starboard  rudder. 
The  process  was  repeated. 

Figure  20  and  Table  7  show  the  results  of  computer 
program  9. 

During  the  spiral  test  of  model  C,  the  angular  velocity 
was  reduced  suddenly  from  r'  =  0.19  to  r'  =  0.28  at  d  =    1.5 
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FIGURE  20.   Spiral  Maneuver  of  Model  C 
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TABLE   VII.       SPIRAL 

MANEUVER  OF   MODEL 

RUDDER   ANGLE 

YAW   RATE 

(DEGREE) 

(RAD./vSEC*L/U) 

:+:  %  ^  +  %  #  s(: 

0.36548072 

-    9.000 

0.35614842 

-    8.000 

0.34629613 

-    7.000 

0.33535161 

-    6.000 

0.32470226 

-    5.000 

0.31269461 

-    4.000 

0.29961270 

-    3.000 

0.28513420 

-    2.000 

o".  26874071 

-    1.000 

0.24949914 

-    0.000 

0.22536594 

1.000 

0.18941998 

2.000 

-0.28294164 

3.000 

-0.29785007 

4.000 

-0.31125194 

5.000 

-0.32350868 

6.000 

-0.33485866 

7.000 

-0.34546858 

8.000 

-0.35545886 

9.000 

-0.36492395 

10.000 

-0.37393415 

11.000 

-0.38254720 

12.000 

-0. 39080832 

11.000 

-0.38356121 

10.000 

-0.37394869 

9.000 

-0.36493850 

8.000 

-0.35547382 

7.000 

-0.34548380 

6.000 

-0.33437481 

5.000 

-0. 32352638 

•    4.000 

-0.31127167 

3.000 

-0.29787230 

2.000 

-0.28296709 

1.000 

-0.26596129 

-    0.000 

-0.24573880 

-    1.000 

-0.21966898 

-    2.000 

-0.17531638 

-    3.000 

0.28510946 

-    4.000 

0.29959112 

-    5.000 

0.31267542 

-    6.000 

0.32468486 

-    7.000 

0.33583581 

-    8.000 

0.34628141 

-    9.000 

0.35613346 

******* 

0.36548072 
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degrees  port  side.   That  means  any  increase  of  the  rudder  angle 
to  port  by  1.5  degrees  will  cause  the  model  C  to  suddenly 
swing  from  starboard  to  port  until  r'  =  0.28.   Likewise,  when 
the  rudder  angle  increase  from  port  to  starboard  the  model  C 
can  turn  against  its  rudder  by  o =   2.5  degrees  starboard  and 
then  suddenly  swing  in  the  opposite  direction  to  a  new  stable 
position  r'  =  0.28. 
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V.   CONCLUSION 

This  thesis  has  presented  the  characteristics  of  the 
dynamically  unstable  ship  in  stability  and  maneuverability . 

As  discussed  in  Part  III,  the  most  powerful  factor  in 
stabilizing  a  unstable  ship  is  a  lifting  surface  installed  aft 
of  midship;  the  addition  of  area  of  fin  is  more  effective 
for  stability  than  the  aspect  ratio.   The  effect  of  a  lifting 
surface  may  be  obtained  by  designing  a  fine  stern  with  neat 
flow  lines  and  deadwood  or  a  proper  rudder.   It  was  found 
that  the  fin  installed  near  the  pivot  point  destabilizes 
the  ship. 

As  discussed  in  part  IV,  the  unstable  ship  is  more 
maneuverable  in  a  turning  maneuver  but  less  maneuverable  in 
the  standard  zig-zag  maneuver  than  the  stable  ship.   In 
turning  maneuver,  the  more  unstable  ship  was  less  sensitive 
to  the  rudder  deflection. 

As  discussed  in  section  IV-B-1 ,  the  followings  are 
considered  necessary  for  optimum  stability  and  maneuverability 
of  a  ship. 

1)  Locating  the  rudder  at  the  point  aft  of  the  midship 
where  the  steady  turning  radius  is  minimum. 

2)  Installing  a  stabilizing  fin  or  deadwood  as  the 
stern  in  order  to  maintain  adequate  stability. 

DSL/360  Language  was  found  to  be  very  powerful  for 
simulating  ship  motion. 
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VI .   RECOMMENDATION 

All  results  in  this  thesis  were  obtained  by  simulation 
using  DSL/360  without  real  model  tests. 

The  author  would  like  to  compare  these  results  from 
real  model  tests  in  the  future.   The  results  in  the  thesis 
are  in  non-dimensional  values  so  that  they  can  be  converted 
to  full  scale  ships  or  models. 

Also  the  following  itmes  are  recommended  for  future  study; 

1)  Maneuverability  of  the  unstable  ship  in  Z-maneuvers 
other  than  the  standard  Z-maneuver. 

2)  Maneuverability  of  a  dynamically  unstable  ship  with 
auto-steering  control. 

3)  Ship  design  for  the  optimum  stability  and  maneuver- 
ability . 
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APPENDIX   A 


Prediction  of  Hydro dynamic  Coefficients  in  Linearized  Equa- 
tion of  the  Ship. 

Using  the  method  in  Ref.  1  and  2,  the  following  predic- 
tions can  be  obtained. 

(1)   A  fin  travelling  at  a  forward  velocity  u  and  trans- 
verse velocity  v  has 


v4  ■ 


3Cl 


ftf  l  3/9    Ji 


(34) 


where 


=  Af/urt 

-    CUjlcl   c& 


Af 
Af 

c, 


(3  =    P/u/f'   o^jAl 


-g 


F 


/+ 


2 


for  general  aspect  ratio  fins 


a 


% 


-   71 


2_         ^rr  for  very  low  aspect  ratio  fins 


(Yr')f 
(A/r')f 


',.//■ 


(34a) 
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(35) 


where 


61     =   aspect  ratio 

^Hjr~   effective  aspect  ratio 

J 


i'f   -  "f/L- 

lis.    =  effective  span 


W*/f 


(Yt)f 


(m 


t 


=  /^y.'f/ 


y;i 


xfJ  ^  /v  /f 


(35a) 


Considering  the  rudder  as  a  fin, 


'i  =  -K 


r 

f 


f       /Y/' 


tf --£<%) 


(35b) 
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(2)   Considering  the  ship's  hull  as  a  fin  with  consider- 
able drag  coefficient  C 


(%')> 


-[(^A+W/ 


(36) 


where 


f^-A -■?«•=" 


cyT\  ='k^+-xf(yj)h 
(hJr)h  =_?»;  ^ '+  x'cyv\ 


£ 


si 


L     .     f 


(36a) 


H   L   z 


hour 


(see   Ref .    1    for  more   details) 
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(3)  Assuming  a  large  deadwood  which  has  a  sufficiently 
low  aspect  ratio,  an  effective  span  equal  to  the  draft  of  the 
ship  and  a  triangle  shape,  than 


'  /V  Jfl  L 

(Yr')fi  =-i(yv')f> 

(^)f,--i-(Yv)fi 

(  /v  \  i  =      4-irAt  hi 


(37) 
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(4)   Assuming  that  a  rudder  or  small  deadwood  has 

"2^      2X 


-n/)      ,   2.      due  to  the  high  aspect  ratio,  then 


2K 


VV*    =-f^v> 


/ 


/ 


M)f>  =i(yv)fz 


<  y*  kz 


47] A {2  hi 


z 


(see  Ref .  1  and  2  for  more  details) 
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(38) 


APPENDIX   B 

Models  for  Stable,  Marginally  Unstable,  Unstable,  very 
Unstable  and  Extremely  Unstable  Ship . 

In  this  section  it  is  intended  to  obtain  models  for 
several  types  of  unstable  ships.   These  will  be  simulated 
using  DSL/360  to  study  their  stability  and  maneuverability 
The  original  model  was  adopted  from  a  five  foot  Series  60 
model  (with  block  coefficient  0.7  and  no  propeller)  as 
originally  developed  by  the  Davidson  Laboratory. 

Principal  characteristics  of  parent  model. 

Length  (L)  5.0  ft. 

Breadth  (B)  0.714  ft. 

Draft  (H)  0.267  ft. 

Block  coefficient  (C  ) 0.7 

LCG  from  bow  (X_) 5.15  ft. 

Displacement  (A) 41. 64  lbs  . 

Area  of  rudder  ( AD) 0  .  021  sq  .  f  t . 

Rudder  span  0.2  f  t . 

Rudder  chord   0 . 105  ft . 

Mass  (m)  1 .  292  slugs 

Mass  coeff.  (m»)  0.2 

Longitudinal  added  mass 

coeff.  (km')  0.004 

Lateral  added  mass  coeff .  (m '  ) 0.18 

2 


76 


Rotational  added  mass  coeff.  (m') 0.165 

2 

Radius  of  gyration  in  air  0.25L 

C.  G.  of  lateral  added  mass  from 

midship  (-£-  )  0  .  024 

Drag  coeff.  at  zero  drift  angle  0.019 


By  using  Appendix  A  and  computer  program  1,  the  parent 
model  can  be  modified  to  the  desired 'models .   From  the 
results  of  computer  program  1,  the  five  models  were  selected, 
which  were  named  stable  model  A,  marginally  unstable  model  B, 
unstable  model  C,  very  unstable  model  D  and  extremely  un- 
stable model  E.   (See  computer  program  1,  Appendix  A  and 
Table  I) 

The  hydrodynamic  derivatives  in  Table  I  were  considered 
to  be  inherent  for  each  model  in  this  thesis. 
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APPENDIX   C 

Selection  of  the  Optimum  Values  of  the  Proportionalities 
of  the  Rudder  Deflection  to  the  Ship's  Heading  error  and 
Yawing  Velocity. 

Let  the  rudder  deflection  equation  be 

Sti)  =  Ai*  +  fa-f  (See   EQ-  23) 

The  ship  motion  under  such  rudder  control  has  the 
characteristic  equation, 


As^Bs^Cs  +  V-O         (seeEq.  25) 


Let 


S   =  'A    T 

n  ,  then 


or 


or       /?%-  pV  b,p  +   h0   -   O  (39) 
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where 


,      Ac  ,      JlZ 

P/~    -fix  and  Uo   ~        ft*  (40) 


From  Eqs.     (24)    and    (25) 

2  =  <r,k+Cz 


(41) 


where 


<J  =  Nii%-»> ')  -//  Mr*7  '** ) 


Likewise 


C    ■=■    CBk,±C+k%+  Cr  (42) 


where 


c. = yv'csJ{-w  %')  -  rfv(Yr'-to') 


c6  i,  («) 
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where 


/  ,  /    .  ' 


cc=(yvtit  -N*Yi) 


Eq.  (41)  gives 

2 


c;-k;+2c,c/^+C 


Eq.  (40)  gives 
Eq.  (40)  gives 


(44) 


8- 


bo  C  C 3 kl  +  3  c*C2  &2%  3  c,  c/  £>2  4-  Cx 3  )  =  /] \  k, 


Thus 


B 


45) 


Inserting   Eq .     (45)    into   Eq.     (44) 


B+fi+Bikl-h&foB,  =o 


where 


33  = 


c )  %b0  /c(A 

3c,xc,.c3  bo 
"c7T 


-  b,  c, 


(46) 
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2. 


In  Eq .  (46)  all  coefficients  except  b   and  b-,  can  be 
obtained  from  hydrod-namic  derivatives  of  a  given  ship. 
Calculator  of  b   and  b-, 

For  3  real  negative  roots  of  the  characterisitic 

3      2 
equator,  As   +  Bs   +  Cs  +  D  =  0 ,  the  point  M  on  Mitrovics 

chart  (see  Ref .  7)  should  be  in  *?  =  1.0  region 


A   (  CuAf  p>&™  ) 


For  3  repeated  roots,  M  should  be  the  cusp  point 
At  the  cusp  point, 


dbo 


=   undefined 


ctbj 

dbo      Jho/Juh 


db,  dbi/ct> 


u)t 


Therefore 


db, 


~-J[      2p-    =  0 
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For  the  third  order  characteristic  Eq.  with   =  1.0 


k  = 


h  :zU)t-2C0tL 


So, 


b,  -  2uJi-3  cOt 


Ah  =  2_  6c0t  =  o 

oia)t 


H      u)±^4; 


b/=    3 


.nd  h0   =   ^7 


(see  Ref.  7  for  more  details) 

Note:   It  is  difficult  to  find  the  exact  cusp  point  by 

the  digital  computer  due  to  the  non-integer 

values  of  b   and  bn  . 
o       1 
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